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I Purpose of Method 

To solve a simultaneous set of n first order differential equations 
y . ■ f (x, y^, yg> ### y n )> given only one initial value for each y. 
atx o . 

II Geometric Principle 

Consider the equation y ■ f (x, y), where (x , y ) is known. Given 

h, find k such that (x + h- y + k) is on the curve. The first 

oo 

approximation to k is given by drawing a line from (x , y ) with 

slope -s^i and finding (x + h, y + k,) on this line. k is 
ox|x«x o o J. £ 

o 

found by picking an intermediate point x- , finding y, on the line, 

then creating a new line through (x , y ) with slope J^i , then 

I x*x, 

finding (x + h, y + k~) on the new line. This process is repeated 
successively until a predetermined number of k. have been found. 
Then k is taken as a weighted average of these k. • 

III Derivation for One Equation 

Given y Q - f (x o> y ), find y such that y « f (x, y). 

7 2>x 3y dx r x I y I 



i it 

y 



2z!l + 2£l $L m f +2ff + ff 2 + (f+ff)f 
3x ^y dx xx xy ""yy v x y J x 3 



The Taylor Series about x is: 

o 



y 



+ < x - X J yl + h < x - x J 2 y" + 7, ( x - x J 3 y"' + ••■ 



o v o' J o 2 1 N o' J o 3 ! 

Let h * x - x Q , 4 y ■ y - y Q > 



(*} A y - f k + i h 2 (f + f f) 1 + i h 3 |f + 2f f + f f 

W u * on 2 x y ' / v 6 L xx xy yy 



( x ~>yJ 



+ (f + f f ) f, 

x y 



n, + ••• 

7 JK> V 
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+ ••• 



Let k x - hf (x o ,y o ) 

k ■ hf(x + mh, y + ink,) 
c o o i. 

k 3 - hf(x o +Ah, y o +^k 2 + (A-(^) kj) 

Where m,X,f are constants. 
Taylor's Formula is: 

*<* ♦ p. y ♦ q) - *<v y ) ♦ f x (* . y ) p + y v y o ) q 

+ 1 [ f xs ( V V P 2 + 2 f *y <V V M + VVo 1 q2 l 
Expanding k 2 , k, by Taylor's Formula, 
k x - h f (x o ,y o ) 

k 2 - h [f (x o ,y o ) + mh f x (x ,y ) ♦ ri^ f y <x ,y ) 

+ I C (mh)2 f xx V.' + ^ k l v w 

k 3 - h {f(x ,y o ) ♦ M x(x .y o ) + [>k 2 + a-p k^fy (x o ,y o ) 

♦ | [ ah) 2 fxx (x Q ,y o ) ♦ 2 A h[/O k2 + U-f) l^yi^) 

By substituting for k., k. and truncating to h , get: 
k x - h f (x o ,y Q ) 
k 2 - h f (x Q ,y o ) ♦ h 2 « (f x ♦ tt) , 

1 (x o»y ) 

+ ± iV(f + 2 f f + f :T) 

2 xx xy yy 



( V y o ) 



hf (V 3 ^ + A ^(f x + y) 



( x ~>yJ 



* I h 3 j^X 2 (* a - + 2 * -* 



XX 



xy 



o"o 



+ f f *) + 2m /> (f + f f ) f 



yy 



"W 
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Form the weighted mean: 

(**) K « ak 1 + bk 2 + ck~ 

where a, b, c are arbitrary constants. 

Forming this sum and setting K ■ A y, ie, equating coefficients in (*) 
and (#*), one gets: 

a + b + c ■ 1 

bra + c^ ■ \ 

bm 2 + c J? - 1/3 

c^>m • l/6 

We now have four equations in six unknowns. 

Upon solving this system for a, b, c, m 9 X * f» get Ay m ak-i + bk« + ck-, 
This agrees with the Taylor expansion of Ay as far as the h term. 
The error is on the order of h . 

IV Ease of Computation 

Four equations in six unknowns gives 2 degrees of freedom. By assigning 
2 values appropriately, one can get particularly simple formulas. A 
widely used set of four k*s is: 

k l " h f ( x o ,y o^ 

k 2 »hf (x + Jh, y + Jkj) 

k 3 - h f (x + £ h, y + \ k 2 ) 
k^ - h f (x + h, y + k^) 

Ay - 1/6 (l^ + 2k 2 + 2k 3 + k^) 

For other simple formulas, see £ if] > P. 186. 

V Generalization to More Equations 

Consider the generalization to n simultaneous differential equations: 
dy i 

d5r " f i (x > y i* y 2> •••» y n } (i a ^ ""> n) 
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Using the simple formula described above, we get: 

l. k u - w 1 (x , yi , y 2 , •••, y n ) 



2. k - M ± (x + | h, y ♦ i k u , y 2 + J k , ...) 

o o 

3. k i3 - hf ± (x Q + J h, y x + i k 12 , y 2 + | k 22 , • •) 

o o 

h. k^ - hf ± (x + h, Tx + k^, y 2 + k 23 , .••) 

o o 

then Ay ± - 1/6 (k ±1 + 2 k ±2 + 2 k ±3 + k^) (i - 1, ••• n) 

VI Modification by S. Gill 

In t2l , S . Gill uses the 2 degrees of freedom to minimize the amount of 
storage required for intermediate computations. If there are n equations, 
Gill succeeds in reducing the order of the number of storage registers 
from Un to 3n (from 8n to 6n for floating point computations). 

VII Accuracy 

To get higher orders of_accuracy, one needs to use a larger number of k's 
in forming the average K. Four k's will give fourth-order accuracy, but 
fifth-order accuracy requires six k's. Every additional k greatly 
increases the number of computations since an additional term in Taylor's 
Formula is required, plus an additional substitution into each f . . 

VIII Recommendation 

This method is particularly useful in obtaining starting values, since 
only one initial condition is required; more sensitive methods can then 
be used to get successive points to greater accuracy. Furthermore, the 
large number of computations makes this method unwieldy for many points. 

The modification by S. Gill is widely used and is recommended. 
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